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invariance are also discussed. 
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1 Introduction 



Nonlinear phenomena in electric conduction play the most important role in many 
electronic applications ranging from single units such as a diode or a transistor 
to entire circuits. For extremely small systems with mesoscopic or atomic length 
scales, such as those which can now be routinely fabricated using nano-technology, 
quantum transport dominates conduction. While we now have a very good under- 
standing of linear quantum transport phenomena in nano-systems where quantum 
coherence plays a vital role, the nonlinear quantum transport properties of meso- 
scopic conductors have received relatively less attention. In this regard, several 
important research results have been reported in recent years[|], |^, ||, H, |^. Ex- 
perimentally, Taboryski et. al.^ have reported observations of nonlinear and 
asymmetric conductance oscillations in quantum point contacts at small bias 
voltages. They found that the non-Ohmic and asymmetric behavior causes a 
rectified DC signal as the response to an applied AC current. On the theoretical 
side several directions have been explored. Wingreen et. a/.0] have presented a 
general formulation to deal with the situation of a non-linear and time-dependent 
current going through a small interacting region where electron energies can be 
changed by time-dependent voltages. De VegvarQ] has studied the low frequency 
second-harmonic transport response of multiprobe mesoscopic conductors using 
a perturbation theory in the framework of Kubo formula and found that the low 
frequency second-harmonic current is a non-Fermi-surface quantity. At the same 
time, Biittiker and his co-workers |1|, |^ have advanced a current conserving the- 
ory for the frequency dependent transport. This theory can be applied to discuss 
the non- linear behavior of mesoscopic samples. It has been recognized H] that 
in non-linear coherent quantum transport, it is essential to consider the internal 
self-consistent potential in order to have the theory satisfy the gauge invariance 
condition. This is a fundamental condition which requires that all physical prop- 
erties predicted by a theory can not change if there is a global voltage shift. 
Recently, Christen and Biittiker [§] have investigated the rectification coefficient 
of a quantum point contact and the non-linear current- voltage characteristic of a 
resonant level in a double barrier structure using this theory of gauge invariant 
non-linear conductance. 

Clearly it is important and useful to further investigate nonlinear quantum 
transport phenomena in coherent quantum conductors. In particular, detailed 
predictions of nonlinear conductance of two-dimensional (2D) systems warrant 
to be made because these systems can now be fabricated in many laboratories. 
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Unfortunately due to various technical difficulties, especially the difficulty of eval- 
uating a quantity called sensitivity (see below), so far the application of Biittiker's 
theory|0 has largely been limited to quasi-lD systems. For 2D conductors, in gen- 
eral some numerical calculation is needed and this has recently been carried out 
by two of the authors to study the admittance of a T-shaped conductor which 
is related to the linear order transmission function. The nonlinear conductance 
for a 2D conductor, on the other hand, has been investigated for a very special 
and exactly solvable model which is a quasi-lD ballistic wire with a (5-function 
impurity confined inside [TT|. Since it is exactly solvable []TD|, the sensitivity 
can be computed in a closed form thereby overcoming the technical difficulties 
associated with the theoretical formalism. To the best of our knowledge, this was 
the only explicit computation of the weakly nonlinear conductance from the gauge 
invariant AC transport theory for a 2D system where mode mixing is the most 
important characteristic. However we note that in order to apply the theoretical 
formalism to a wide range of 2D mesoscopic conductors, a more general numeri- 
cal method must be developed and various physical issues cleared. The purpose 
of this article is to report our development of such a numerical method, and to 
investigate the weakly nonlinear transport properties of a truly 2D conductor. 

As we have noted from the previous investigation of the exactly solvable 
model|T^, for a geometrically symmetric system the second order non-linear con- 
ductance Gill niust be zero from a general argument (see below). Hence the 
non-linear effect, i.e. a non-zero Ga/s^, obtained in Ref. []10| is a delicate ef- 
fect of the asymmetric scattering boundary||12||. Such an asymmetry is brought 
about when the (5-function scatterer is not located at the center of the scattering 
volume [|1^]. Already, very interesting and physically revealing behavior of the 
local current response (the sensitivity) has been found. In this work, on the other 
hand, we shall focus on a much more general situation by investigating the 2D 
conductor depicted in Fig. (1) where the scattering volume is defined by the 
shaded area. The two leads with width W are assumed to extend far away from 
the scattering volume. The shape of the side stub is controlled by the parameter 
H as shown in Fig. (1), hence various different 2D systems can be generated by 
varying H. For H = W, the scattering volume is a geometrically asymmetric sys- 
tem, but the asymmetry is only due to the asymmetric locations of the scattering 
volume boundary. For this case we thus expect that the physics should be similar 
to those obtained in Ref. |10|. For H = 2W, the scattering volume becomes a 
geometrically symmetric system where Ga/s^ must be zero. For other values of 
H between W and 2W, the scattering volume is intrinsically asymmetric. By 
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varying H, we shall study the crossover behavior of Gm between the symmetric 
and asymmetric situations. 

Our results show that the external (the internal) response of the second order 
non-linear conductance changes sign from negative (positive) to positive (neg- 
ative) near a quantum resonant point. The cancellation of the external and 
internal responses results in a much smaller second order non-linear conductance 
Gill, i-G., Gill is one order of magnitude smaller than the external or internal 
response. The behavior of Gm is non-monotonic when changing the parameter 
H in the range W < H < 2W: this is because Gm is very small at H = W 
as it is solely due to the asymmetric scattering boundary, it increases as H is 
increased, and it is zero aX H = 2W. Another result of our analysis concerns 
the gauge invariant condition J2j Gap-y = 0. It turns out that for systems with 
a finite scattering volume as those of any numerical calculations, if the global 
partial density of states (GPDOS) is computed from the energy derivatives of 
the scattering matrix, it was found ||10| that a correction term must be added to 
satisfy the gauge invariant condition. For the exactly solvable model studied in 
our previous work[|10[, this correction term has been derived[0] analytically. We 
shall examine this effect for the conductor studied here. 

The paper is organized as follows. In the next section we shall briefly review 
the gauge invariant theory for non-linear transport developed by Biittiker[0. The 
method of calculating various quantities needed for non-linear conductance and 
our results are presented in sections 3 and 4. The last section serves as a brief 
summary. 



2 The formalism 

The gauge invariant formalism of nonlinear transport has been developed and 
clearly discussed in Ref. and we refer interested reader to the original work. 
In this section, we shall outline the main steps of the application of this formalism 
for our calculation. For a multi-probe mesoscopic conductor, the current through 
probe a is given by|l|, |^ 

= X E / dEf{E -Ef- eVp)A^piE, {V,}), (1) 

where f{E) is the Fermi distribution function, and 

A^^iE, {V,}) = Tr[lJ^^ - sipiE, {V,})s^piE, {V,})] (2) 
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are the screened (negative) transmission functions which are expressed in terms 
of the scattering matrix s,^^. For weakly non-hnear transport, Eq.(|I|) can be 
expanded with respect to the voltages V^. For a conductor which has only two 
probes, up to the second order nonlinear term such an expansion leads to the 



following equation p, |10 



h = Gn{V,-V2) + Gni{V,-V2y , (3) 

where Gu is the usual linear conductance and Gm is the second order nonlinear 
conductance which we wish to compute for the conductor of Fig. (1). 

It can further be proven that Gap-y is the sum of two terms which are the 
external and internal contributions: 

Gap-i = G^^^ + (4) 

where the external contribution can be obtained using the free electron scattering 
theory by evaluating the energy derivatives of the scattering matrix, 

2 

Glp^ = jI dE{-dEf)edEA^p5p^ . (5) 

The internal contribution, on the other hand, is much more difficult to obtain 
because it depends on the potential derivatives of the scattering matrix, 

2 

Gif3y = jldE{-dEf)idv,A^p + dv,A^^) . (6) 

The reason why this is difficult to evaluate is because when the voltage of a probe 
changes, the entire potential landscape of the scattering volume will change 
accordingly through the electron-electron interactions. Hence the internal con- 
tribution to the nonlinear conductance can be obtained only after an interacting 
electron problem has been solved]^. This is a very difficult task and so far has 
not been successfully implemented in a numerical scheme. However if we can use 
the Thomas-Fermi linear screening model, which is more appropriate for metallic 
conductors, the internal contribution can be computed through the evaluation 
of quantities called sensitivity and characteristic potential^]. It can be shownp 
that the potential derivative of the transmission function is given as, 

dv^A^p = ^T^ J dhr]ap{r)u^ir) (7) 

where 



1 SA^p ^ 1_ I 6sap 5slf^ . 

47r 5U{r) ATI "^^^"^ 5t/(r) ^ 5f/(r) ' 
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is called sensitivity^^ which measures the local electric current response to 
an external perturbation, ^-^(r) is the characteristic potential which measures 
the variation of the potential landscape of the scattering volume due to the 
perturbation |T]]. Within the Thomas- Fermi screening model, it is given by 

_ dn{r,j) dnjr) 
'~ dE ' dE ■ ^ ' 

Here the local partial density of states (LPDOS) dn{r,'-f)/dE is called the injec- 
tivity and is given by|T^ the scattering wavef unctions, 



dE ^ hv^ ' ^ ^ 



where v^n is the electron velocity for the propagating channel labeled by n, and 
dn{r)/dE = J2a dn{r, a)/dE is the total local density of states. 

From the weakly nonlinear conductance formalism summarized above, several 
observations are in order. First, from an application of this formalism point of 
view, a crucial step is the evaluation of the sensitivity rjap which depends on the 
functional derivative of the scattering matrix with respect to the local potential 
variation. The latter is caused by the external perturbation, i.e. the change of 
the electrochemical potential at a lead. We are aware of two ways of calculating 



the sensitivity [T^]. The first is to evaluate dSap/SU directly by introducing a 6- 
function of infinitesimal strength 6U inside the scattering region. Alternatively, 
one can calculate it using the retarded Green's function. For a 2D system, in 
general the Green's function can not be obtained explicitly except in very special 
cases such as that studied in Ref. |ll|], hence we shall use the first method by 



directly computing the sensitivity. As a second observation which is physically 
important 1 112|| , we can discuss the general behavior of the nonlinear conductance 
Gill. From Eq. (^, for a symmetric scattering volume with scattering potential 
U{x,y) = U{—x,y) where x is the propagation direction, we must have — Ji if 
Vi and V2 are interchanged. Hence we conclude that for a symmetric scattering 
volume there is no quadratic terms, i.e., Gm = 0. On the other hand, in general 
Gill 7^ for geometrically asymmetrical systems. Finally, due to the current 
conservation and the gauge invariance condition, namely the entire physics is 
independent of a global voltage shift, it is not difficult to prove ||, [l^ 



— X! — X! ^Q/37 — . (11) 

a /3 7 



Our results will allow a direct confirmation of this equation. 
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3 Numerical method 



There are several ways to solve the scattering matrix of 2D ballistic conductors, 
such as the mode-matching method[|l^, the recursive Green's function method ||T6|, 
and the finite-element method |I9|, |20|. However we found that all these 



methods are not particularly easy to apply here because we need not only the 
scattering matrix, but also the sensitivity rjap- For this purpose, we found that a 
technique for computing scattering matrix which is developed in Ref. ||2l| is quite 



useful and we shall discuss, in some detail, our numerical procedure for finding 

In particular, we construct a global scattering matrix using the mode-matching 
method of Ref. ||2l|. For a scattering volume which is not uniform along its 
longitudinal direction, we divide it into a number of uniform sections, e.g., the 
asymmetric cavity as shown in Fig. (1) can be divided into four uniform sections. 
The scattering matrix associated with the n-th section S„ is the composition of 
two individual scattering matrices and S^, i.e., S„ = S{ ® where <S) is the 
operator which denotes the composition of two scattering matricesp2|. Here S{ 



describes the free propagation from the left end of the n-th section to its right 
end and is given by 

P, 







(12) 



gjfc™L„^ jg ^Yie longi- 



where Pn is a diagonal matrix with elements (Pn)mm, 
tudinal wave number for the m-th mode and L„ is the length of the n-th section. 
The scattering process at the interface between two adjacent sections (the n-th 
and (n + l)-th sections) is described by S^. Care must be taken when matching 
the wavefunctions of two sections with different widths at the section boundary. 
If the width of the n-th section Wn is not greater than Wn+i, we have 



I 

CK 



n+1 



I 



CK 



n+1 



(13) 



where K„ is a diagonal matrix with diagonal element k^, and I is a unit matrix. 



C is a matrix which denotes the coupling between the transverse modes in the 
two sections and its elements are given by Cij = {(pn\(pn+i) where 0*„ is the ith 
transverse mode in the nth section. C'^ is the transpose of the matrix C. On the 
other hand if Wn > Wn+i, we have 



-I 

C^K. 



C 



I 

C^K. 



-C 

K„,+i 



(14) 
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Once the scattering matrix for each section is known, the global scattering ma- 
trix can be easily constructed by the composition of all the individual scattering 
matrices, 

S = Si ® S2 ® . . . ® Sm-1 (15) 

where M is the total number of sections. 

It should be noted that the global scattering matrix S calculated this way 



is different from the standard one|T^ which only connects the outgoing wave 
to the incoming wave. Here, the global scattering matrix S involves the non- 
propagating channels as well and does not satisfy the unitarity condition. In 
order to obtain a physical scattering matrix which we denote by the lower case s, 
we first rewrite global S in the form of 2x2 subblocks and obtain four sub-matrices 
Sjj where = 1,2). Then for each sub-matrix Sjj we build a new matrix s^- 
constructed by the first Nq rows and columns of Sjj, where Nq is the number 
of propagating channels. By writing the four newly constructed matrices in the 
form of 2 X 2 blocks, a 2Ai"o-dimensional scattering matrix s is obtained which 
is the true scattering matrix that connects the outgoing wave to the incoming 
wave. In order to obtain a unitary scattering matrix s, one should further take a 
unitary transformation. 



AsA-i 
V 
V 



(16) 



where V is a A^^q- dimensional diagonal matrix with diagonal element yk^. 

The above procedures can easily be modified to compute the sensitivity rjapijo). 
For this purpose, we shall make use of a 5-function impurity to calculate the 
functional derivatives of the scattering matrices dsap / 5U {t) . This is achieved 
as follows. We put a 5-function impurity with infinitesimal strength 7, V^(r) = 
7(5(r — To), at arbitrary positions r = tq in the scattering volume. We then 
calculate the scattering matrix as a function of 7. Finally we use a five-point 
numerical derivative to evaluate dsap/ 5U{r:) = dSa/3/d-y\j=o. With this result we 
can obtain the sensitivity from Eq. (^. 

The scattering matrix can still be obtained using the approach discussed above 
even including the ^-function impurity. In fact it has been derived in the presence 
of a magnetic field by Tamura and Ando||22|. Here we give the expression in the 
absence of the field[^. Suppose the 5-function impurity is located in the n-th 
section at position fq = (xq, yo), where xq and yo is the distance from the left and 
bottom boundary of the section respectively. The scattering matrix associated 
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with this section is then given by 



S„ = S{(xo) ® ® S{(L„ - xo) ® . 



(17) 



Here describes the scattering process associated with the 5-function impurity 
and is given by[p3| 



-I I 



I -I 

iKn + r iKn 



(18) 



where the matrix T describes the mode-mixing effect due to the 5-function impu- 
rity and its matrix elements are given by Tpg = 27 sm{pTTyo/Wn) sin(g7r7/o/PVn) /Wn 
with Wn being the width of the section. With the ^-function included this way, 
we can again apply our method, Eq. (|15D , to compute the scattering matrix 
s = 3(7) and complete the numerical derivatives discussed in the last paragraph. 

To end this section, we briefly mention two other points. First, the character- 
istic potential as given by Eq. (P) is evaluated using the scattering wavefunction 
which can be calculated in two ways, directly or indirectly. One can directly 
compute the wavefunction using the mode matching method ||15[| or flnite element 
method ||T9|, pO| . Or one can compute the wavefunction (|^l/p) indirectly by com- 
puting the local partial DOS called emissivity defined as[jl4i 



dn{a, r) 
dE 



1 



Ss 



0/3 



47rz ^" " ^ "^5ef/(r; 



- s, 



a/3 



6eU{r) 



(19) 



The microreversibility of the scattering matrix implies, 

dn{a,r,B) dn{r,a,—B) 



where B is the magnetic field. We will use Eqs.(|lO]),(|l9D, and (|20|) to compute the 
wavefunction |\Efp since we have to compute the scattering matrix or 5s/5U{y) 
anyway for the sensitivity. Second, the energy derivatives of the transmission 
function, which determines the external response contribution to the second order 
nonlinear conductance, is evaluated using a five point numerical difference. These 
procedures are the same as that of our earlier work[|]. 



4 Results 



As a first result we plot in Fig. (2) the sensitivity 7711 (r) at three different po- 
sitions inside the scattering volume, ri = (0.5W, 0.5iy), Y2 = {0.5W,1.25W), 
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and rs = {1.25W,0.625W), as a function of the normalized electron momentum 
kpW/TT where kp is the electron Fermi wave number. These results are for an 
asymmetric system (see Fig. (1)) with the parameter H = 1.2514^. As discussed 
in Section 2 rja/B appears naturally in the theoretical formalism, and it essentially 
describes the local (internal) electric current response of the scattering problem 
when there is a small local potential change. It is related to the real part of the 
diagonal elements of the scattering Green's function ||13||. From Fig. (2), we see 
that different positions inside the scattering volume have quite different inter- 
nal responses in terms of the peak heights of the apparent resonance behavior. 
On the other hand, the peak positions occur at the same electron energies given 
by kpW/Ti for rju at all three positions. We have checked (see below and Fig. 
(5b)) that the peak positions also coincide with those of the conductance. Hence 
we may conclude that the local current response can have sharp changes, from 
positive values to negative values, across the energy of a resonance which also 
mediates a resonance transmission. For the two positions ri and r2 which are 
located in the left part of the cavity, the shapes of the sensitivities are more sim- 
ilar to each other. This is to be compared with that of the position which is 
located in the right part of the cavity. The differences are evident from the three 
curves of Fig. (2). 

To get a more intuitive picture of the spatial dependence of the sensitivity, 
in Fig. (3) we plot ?7ii(r) in the entire scattering volume for two different values 
of the electron momentum {H = 1.25W). First is for kplV/TT = 1.715 which 
is off resonance, while the second is for kplV/ir = 1.795 which is on resonance. 
From the lower panel of Fig. (3) which corresponds to the resonant energy, the 
behavior of rju is reminiscent of a standing wave which is in accordance to our 
usual picture of a quantum resonance. The positions ri and r2 of Fig. (2) are 
located at a peak of the sensitivity profile while is at a valley. This explains 
why in Fig. (2) we observe the large resonant peak at ri and r2 but not at 
r3. When off resonance, the upper panel of Fig. (3) shows less regular patterns 
for rjii. Hence the local current response can behave regularly or less regularly 
depending on the electron Fermi energy being on or off a quantum resonance of 
the scattering cavity. In comparison, for both ID and 2D scattering problems 
involving a 6 potential barrier, the sensitivity has been derived analytically in 
Ref. and |jlO[. There, rjn shows strong spatial regular oscillations. 

When the geometry parameter H = W, the conductor becomes a T-shaped 
junction. The upper panel of Fig. (4) shows rju for this situation at kplV/n = 
1.325. The lower panel of Fig. (4) plots the characteristic potential 2ui{r) — 1 for 
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this case. The quantity 2ui{r) — 1 is interesting because it can easily be shown[||, 
using Eqs. (|^), (|^), (|^) and applying the gauge invariance condition (0), that 
the nonlinear conductance can be re-written as 

Gai3y = 'iTr^J dE{-dEf) J dh{r]apu^{r)+riayUi3{r)-riaf3S^p) . (21) 

Hence the quantity 2ui{r) — 1 appears naturally in this form of Cm. From Fig. 
(4) it is clear that rju is symmetric and 2mi — 1 is anti-symmetric along x axis. As a 
result, Gill will be zero for this T-shaped junction if the scattering volume is sym- 



metric due to the spatial integration of Eq. (pi]). To systematically investigate 
the behavior of Gm as the conductor shape changes from symmetric to asymmet- 
ric, we have calculated this quantity for several values of the geometric parameter 
H at zero temperature: H = W,1.25W,1.5W , and 1.751^ (see Fig. (1)). Figs. 
(5a-5d) plot the DC-conductance On, the external and internal responses of the 
second order non-linear conductance, G\^i, and Gm as a function of the 

normalized electron momentum kpW/Tr for these configurations. Several interest- 
ing features have been observed. First of all, the external (the internal) response 
of the second order non-linear conductance changes sign from negative (positive) 
to positive (negative) near the resonant point. This behavior is similar to that 
of one-dimensional asymmetric double barrier resonant tunneling In that case 
Gill = {e^/h){dT/dE){T2 — ri)/r, where T is the transmission coefficient, Fj is 
the decay width of each barrier, and F = Fi + F2. Because of the term dT/dE, 
Gill changes sign across the resonant point and hence can be negative. The can- 
cellation of the external and internal responses results in a much smaller Gm: 
one order of magnitude smaller than the internal or external contribution alone. 
Secondly, for H = W the asymmetry of the scattering volume only comes from 
the location of the scattering volume boundary which is at our disposal, Gm has 
the smallest values for all H < 2W studied. For W < H < 2W, the conductor 
is intrinsically asymmetric and Gm are larger. Also the resonance behavior of 
Gm becomes substantially sharper as H is increased. While Gm increases as H 
increases from H = W, it eventually starts to decrease after H > 1.51^. This 
is clearly seen from Fig. (5d) for which H = 1.75W. This is because for larger 
H the system approaches to the symmetric conductor at H = 2W, for which 
Gm = as discussed above. 

So far, as noted in Ref. we have computed the external contribution 

Gli-f^ of Eq. d^) using a procedure which employs the gauge invariance condition. 
However if we directly use the right hand side of Eq. (H) to compute Gl^, the 
result will not be accurate because the scattering volume is finite. In particular, as 



11 



studied previouslyfT^ for a symmetric system Gm will be non-zero if calculated 
using the right hand side of Eq.(^. In terms of the gauge invariant condition, 
this will lead to a violation, i.e. Gm + Gu2 7^ 0. While such an error is not 
relevant if the scattering volume is very large, for our numerical calculations it 
will generate incorrect conclusions since the scattering volume is always finite. 
Hence, in order to use the right hand side to compute Gin, ^ correction term is 
needed to preserve the gauge invariance. 

For a quasi-lD wire with a 6 function impurity, this correction term has been 
derived|T^ analytically. For that situation, the correction, C, consists of two 
terms ||lO|: 



C= ^/?e(sn) + i?e(^^iJ^e*^"("^-^^)) . (22) 
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where fc„ is the longitudinal momentum for the n-th mode with k"^ = kp — 
{riTc/Wy and E = h'^kp/2m, Xi and X2 are the coordinates for the scattering 
boundaries, and 6„ is the reflection scattering amplitude. Clearly, the first term 
is oscillating as the linear size of the scattering volume increases (due to Sn) and 
it is only relevant near the edge of the first propagation threshold; the second term 
is exponentially decaying to zero as the size of the volume increases and it comes 
solely from mode mixing and contributed by the evanescent modes. Although this 
form of the correction term was derived from another system, it is nevertheless 
interesting to compare this formula for the conductor studied here. 

In Fig. (6) we plot the correction term Eq. (|22|) together with Gm + G112 
which were computed using the right hand side of Eq. (^, for the case H = W. 
Eq. (^) was evaluated using the scattering matrix elements obtained from our 
numerical calculation, and the system size {x2 — xi) was specified as the length of 
our scattering volume 2L + 2W (see Fig. (1)). For this conductor with H = W, 
we expect a good comparison with Eq. (p^) which was derived for a 5-function 
scatterer inside a wire, because for both systems the geometric asymmetry is 
solely due to the positions of the scattering volume boundary. Other than that, 
these systems are actually symmetric with respect to the scattering potential. 
Indeed, Fig. (6) clearly shows that there is essentially no difference between Eq. 
( p2D and our numerical data when the scattering volume is large (L = 4W). 
On the other hand, for a conductor with H = 1.5W which has an intrinsically 
asymmetric scattering volume, the comparison is qualitative as shown in Fig. (7). 
However the trend of the two curves are still similar. We may thus conclude that 
for the gauge invariant condition, the correction term to the external contribution 
Eq. d^) has a form with the same nature as that of Eq. (|2^ ) above. 
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5 Summary 



In this work we have developed a numerical technique based on a scattering ma- 
trix to compute the weakly nonlinear conductance. This technique is particularly 
useful for conductors whose scattering volume can be naturally divided into sev- 
eral regions. The most difficult step is the evaluation of the local electric current 
response, namely the sensitivity. We have reported how to obtain this quantity 
numerically, thus further investigations of the interesting nonlinear conductance 
problem can be carried out using our numerical method. We have found that 
sensitivity behaves differently when the transport energy is on or off resonance. 
The former leads to standing-wave type spatial dependence, while the latter is 
behaving in a less regular fashion. In all cases, the sensitivity shows spatial os- 
cillating pattern, which is similar to those known from exact calculations for ID 
models. 

The nonlinear conductance can be non-zero only for geometrically asymmetric 
systems. The asymmetry can be introduced in two ways. The first is through 
the intrinsically asymmetrical shape of a conductor such as those of Fig. (1) 
with W < H < 2W. The other, which is a trivial asymmetry, is through the 
asymmetrical location of the scattering volume boundary, e.g. the case H = 
W. We discover that the intrinsical asymmetry leads to much larger nonlinear 
conductance than the other case. Furthermore, for the symmetrical scattering 
junction but with asymmetrical location of the boundary {H = W case), for large 
size L the behavior of the gauge invariance condition agrees almost perfectly with 
Eq. which was derived from a completely different system but also with 

the asymmetry introduced by the location of the scattering volume boundary 
only. On the other hand, such an agreement is less perfect for the intrinsically 
asymmetrical system. Hence we may conclude that the nonlinear conductance 
behaves in quite different manner depending on how the asymmetry is introduced. 

The sign of the nonlinear conductance can be positive or negative. Very sharp 
variations of this quantity is discovered at quantum resonances for the conductor 
studied here, where such resonances are marked by sharp reductions of the linear 
conductance Gu. Hence near a resonance, the electric current may actually 
decrease for an increasing voltage difference by Eq. since Gm is negative. 
Such a behavior is precisely the expected nonlinear conduction characteristic, and 
up to the second order in voltage difference our results can provide a prediction. 
Clearly, as the voltage difference becomes large, even higher order conductances 
must be included in order to have a meaningful prediction of the nonlinear I-V 
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curve. 
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Figure Captions 



Fig. 1 Schematic view of an asymmetric cavity (the shaded area) embedded in a 
quantum wire. 

Fig. 2 The sensitivity rju at three different positions, (x, y)={(0.5W, 0.5W), 
(0.5W, 1.25W), (1.25, 0.625W)}, as a function of the normahzed electron 
momentum kFW/ir for H = 1.25W. 

Fig. 3 Three dimensional view of the sensitivity rjn for H = 1.25W at two different 
values of the electron momentum kplV/TT = 1.715 and kpW/TT = 1.795. 

Fig. 4 Three dimensional view of the sensitivity rjn and the characteristic potential 
2ui — 1 for a symmetric T-shaped cavity at kpW/TT = 1.325 and H = W. 

Fig. 5 The DC-conductances Gu and the leading order nonlinear terms Gm as a 
function of the normalized electron momentum kpyV/TT, solid lines for Gin 
and dotted lines for Gl^. (a) H = W, (h) H = 1.25W, (c) H = 1.5W, 
(d) H = 1.75W. Here Ei is the threshold of the first subband defined as 
E, = h\y{2mW^). 

Fig. 6 Gill + 6*112 (solid line) and the correction term according to Eq.(p2D (dotted 
line) versus momentum for the T-shaped cavity {H = W). Upper panel: 
L = 0. Lower panel: L = 4W. 

Fig. 7 Gill + 6*112 (solid line) and the correction term according to Eq.(^2|) (dotted 
line) versus momentum for an asymmetric structure {H = 1.5W). Upper 
panel: L = 0. Lower panel: L = 4W. 
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